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ABSTRACT 

Wc have used multi-band imaging to investigate the nature of an extreme starburst 
environment in the nearby Lyman-break galaxy analogue Haro 11 (ESO 350-IG038) by 
means of its stellar cluster population. The central starburst region has been observed 
in 8 different high resolution T/S'T wavebands, sampling the stellar and gas components 
from UV to NIR. Photometric imaging of the galaxy was also carried out at 2.16 //m 
r ) ' by NaCo AO instrument at the ESO Very Large Telescope (VLT). We constructed 

^*^ I integrated spectral energy distributions (SEDs) for about 200 star clusters located 

r^ • in the active star-forming regions and compared them with single stellar population 

Mh' models (suitable for physical properties of very young cluster population) in order to 

Q , derive ages, masses and extinctions of the star clusters. The cluster age distribution we 

$— ( ■ recover confirms that the present starburst has lasted for 40 Myr, and shows a peak of 

tj^ I cluster formation only 3.5 Myr old. With such an extremely young cluster population, 

C3 . Haro 11 represents a unique opportunity to investigate the youngest phase of the 

cluster formation process and evolution in starburst systems. We looked for possible 
^vj ! relations between cluster ages, extinctions, and masses. Extinction tends to diminish as 

^ . function of the cluster age, but the spread is large and reaches the highest dispersion 

QQ ' for clusters in partial embedded phases (< 5 Myr). A fraction of low-mass (below 

l/^ i IO^Mq), very young (1 — 3 Myr) clusters is missing, either because they are embedded 

\^ • in the parental molecular cloud and heavily extinguished, or because of blending with 

neighbouring clusters. The range of the cluster masses is wide, we observe that more 
ly-^ . than 30 % of the clusters have masses above 10^ Mq, qualifying them as super star 

f~^ ■ clusters. Almost half of the cluster sample is affected by flux excesses at wavelengths 

f^ I > 8000 A which cannot be explained by simple stellar evolutionary models. Fitting 

SED models over all wavebands leads to systematic overestimates of cluster ages and 
incorrect masses for the stellar population supplying the light in these clusters. We 
show that the red excess affects also the HST F814W filter, which is typically used 
to constrain cluster physical properties. The clusters which show the red excess are 
^ ' younger than 40 Myr; we discuss possible physical explanations for the phenomenon. 

Finally, we estimate that Haro 11 has produced bound clusters at a rate almost a 
factor of 10 higher than the massive and regular spirals, like the Milky Way. The 
present cluster formation efficiency is ~ 38% of the galactic star formation rate. 

Key words: galaxies: starburst - galaxies: star clusters: individual: Haro 11, - galax- 
ies: irregular 



1 INTRODUCTION SCs are still, or in most of the cases, embedded in their 

dusty cocoons and at constant risk of destruction by in- 
Youns; star clusters (SCs) are luminous and quite compact , ,. .,, ,i . , , i, i- ■ , i r n i 

^ ' ,,,, ,, teractions with the interstellar medium, inte rnal leedback, 

obiects which m some ways resemble older massive elob- ^- ^ ^ r ^ i- / [T ||r,r,i A j r 

, , , , ,f tidal lorces, and evaporation (see lLarsenll201CI and reierences 

ular clusters. However while these old systems are pas- ,, . \ m, , , i ,, i , ■ ,.,i 

therein), the youngest phases, when the cluster is still par- 
sively radma' due to stellar and secular evolution, youne , . ,, i , i , ■ ., , , i i , i ^-m ,. 

tially embedded m its natal molecular cloud, are dithcult 
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to explore because of the high extinction level s which intro- 
duce uncertainties in age and mass estimates. iLada fc Ladal 
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l|2003l ) observed that the embedded phase in young Milky 
Way (MW) clusters last s a few million years. Ver y young 
clusters in the Antennae (jWhitmore fc Zhaneil2002f ) system, 
with ages ^ 4 Myr, sho w a wide range of extinctio n between 
0.5 ^ Av ^ 7.6 mag. IWhitmore fc Zhand [20o3 identified 
an additional population of "red clusters", which they in- 
terpreted as young embedded clusters. Among them, the 2 
Myr old WS 80, with its extreme V — I = 2.92 colour, mass 
4 X l(f Mq and Ay = 7.6 ma g, is the brightest infrared 
source in the Antennae system. IWhitmore fc Zhang (12003) 
indicated roughly 6 Myr as the age at which a young clus- 
ter has gotten rid of the dusty material left over from its 
parent cloud. In another galaxy, NGC 4449, iReines et al.l 
l|2008al ) found 3-5 Myr old clusters with relatively low ex- 
tinctions {Av ^ 1.5 mag). Although these very early stages 
are fundamental for our understanding of cluster formation, 
evolution and possible destruction (90 % of cl usters are not 
expe cted to survive after their first 10 Myr, iLada fc Ladal 
120031 ). they are still poorly understood. 

The cluster formation process is intimately connected 
both to the galactic environment via the star formation 
efficiency (SFE), and on long time scales, to the global 
star formation history (SFH) of the galaxy. Interactions 
and/or merging between gas-rich galaxies provide physi- 
cal conditions which favour the formation of star clusters 
and enhance star formation rates (SFR). Star clusters are 
easily detected in active star fo r ming regions inside nor- 
mal spiral galaxies (jLarsenl |2002|: ICantiello et al.ll2009l): in 
starburst galaxies, like M82 (Konstantopoulos ct al. 20091 
and references therein ) and NGC 3310 (dc Griis ct al. 2003a,; 
Ide Griis et al.ll2003bl); interacting jBastian et al.ll2005l') and 



merg i ng systems ( Whitmore et al.l 1 19991 : IVaisanen et al.l 
120081 : IWilson et al.l 120061 ) to cite only a few examples. In- 



terestingly, very bright and massive clusters have also been 
detected in smaller and less regula r systems, like the d warf 
post-starbur st galaxies NGC 15 6 9 (lAnders et al. J200J) and 
NGC 1709 (|Billett et al.1 |2002| : lAnnibali et all I2OO9I ) and 
in blue compact galaxies (BCGs), whic h are active dwarf 
starburst galaxies with low metallicity (jOstlin et al.ll2003l . 
lAdamo et al.ll2009l : Adamo et al. in prep.). 

The cluster initial mass function (CIMF) and the cluster 
luminosity function (CLF) are directly connected to the star 
formation processes which act from stellar to galactic scales, 
due t o the clustered nature of star formation ( Lada fc Ladal 
I2OO3I ). In general, the CIMF can be approximated by a 
power la w relation of the type dN(M)/dM = CM~~' , where 
7 » 2 j Zhang fc Falll 1 19991 : iBik et al.l l2003l : iHunter et all 
|2003| . etc.). The CLF, dN{L) oc L'^'dL, on the other hand, 
has been shown to have wider range of the index values, 
1.8 ^ a ^ 2.7, and in some cases is better fit with a double 
power law with a break in the range —10.4 ^ My ^ —8 mag 
(|Larsenll2003 : IWhitmore et al.lll999l : lGieles et al1l2006bl V In 
a galaxy, the CLF is the synthesis of several populations of 
SCs assembled at different ages. For this reason, the CLF 
contains imprints of both the mechanisms behind cluster 
evolution and the initial environmental conditions where the 
clusters form. A number of studies have shown that high 
star formation levels enable a better sampling of the LF 
up to its brighter end and in turn d etermine the formation 
of the brigh test cluster in a galaxy (|Larsen ll200J: Gieles et 
al. 2006a.b: [BastianI l2008l : iLarsenl [20091 : ICiclcs 200i). On 
the other hand, physical factors also put constraints on the 



mass (and consequently the luminosity) of the clusters which 
are born, such as mass of the parent giant molecular cloud 
and the pr essure a nd density of the interstellar medium 
jLada fc Lada^OOJ : lBillett et a"Lll2002l . iHunter et aLll2003l '). 
ICieles et al.l (|2006al ) noted that if physical conditions in the 
galactic environment limit the maximum value of the cluster 
mass, leading to a truncated CIMF, this effect will appear 
in the LF as a break which corresponds to the maximum 
cluster mass at the oldest age the galaxy can form. 

In this paper we focus our attention on the blue com- 
pact galaxy Haro 11 (ESO350-IG038), which has been 
shown to be extremely efficient in producing SCs. Gen- 
erally, BCGs have magnitude betw een Mb ~ —12 and 
Ms « -21 JKunth fc Ostlinl I2OO0I ). They are character- 
ized by an ongoing active starburst phase with high star 
formation rates (SFRs), producing their typical high sur- 
face brightness and prominent nebular emission line spec- 
tra. Photometric studies of their galactic stellar popula- 
tions have reveal ed the presence of a n old underlining stel- 
lar componerit jKunth et al.1 1 19881 : iBergvall 6s"tiinl I2OO2I : 
iMicheva et al.ll2010l ). Thanks to the HSVs resolution power, 
it has been possible to resolve the BCG burst regions into 
many bright star clusters implying that the star formation is 
remarkably efficient. However, it still not clear yet how sys- 
tems with mass es g^ 10^° Mp) can f orm such large numbers 
of star clusters. iBillett et al.l (|2002h . investigating a sample 
of 22 nearby (< 7 Mpc) dwarf starburst and post-starburst 
galaxies, noticed that dwarf starburst galaxies have lumi- 
nosities brighter than My = — 16 and pointed out that de- 
spite their sizes, they are able to host very massive clusters. 
Whether this is a consequence of stochastic sampling effects 
of the CIMF in combination with environmental conditions 
which favour the formation of star clusters with M ^ 10^ Mq 
or different cluster formation physics is not yet established. 
In this paper, we will call clusters with masses ^ W^Mq 
super star clusters (SSCs). 

In Figure |TJ we show a three-colour image of Haro 11. 
The t hree main st a rburst regions, called knots A, B, and 
C by iKunth et al.l (|2003l l , are easily traced in the heart- 
shape. Dust filaments cross knot B and are present in all 
the burs t regions. Haro 11 is a Lya and Lyman continuum 
emitter dHaves et al.l 20071) a nd can be considered a local 
analog (|Overzier et al.l [2003) of the high redshift Lyman 
break galaxies (LBGs). Moreover, Haro ll's IR luminosity 
of 1.9 X 10^^ L0 (section 6.2) locates Haro 11 in the rang e 
of the luminous IR galaxies (LIRGs). lOstlin et al] (|2001r l. 
studying the Hq velocity field, found that the mass esti- 
mate from rotation curves is lower than the photometric 
one, probably because the galaxy is not rotationally sup- 
ported, perhaps because the gas is not in dynamical equilib- 
rium. The total estimated stellar mass is ~ IO^^Mq. Haro 
ll's perturbed morphology, multiple-component Ha veloc- 
ity field and high SFR are all signatures of a merger be- 
tween a low-mass, evolved system an d a gas-rich compo- 
nent. Interestingly. lOstlin et al.l (|200H ) have estimated that 
the present burst in Haro 11 has been active for ~ 35 Myr, 
and using the mass of neutral hydrogen. Hi estimated by 
IBergvall et al.l ((2000|), a gas consumption time scale (that 
is, remaining lifetime) of ~ 5 Myr for the present burst 
phase. We estimated in section 16.21 that the present SFR 
in Haro 11 is of 22 ± 3MQyr~^ . If we take into account the 
total mass of gas in Hi, Hll, H2 and photodissociation re- 
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gions, M]j°* ~ 2 X IO'^Mq (JBergvall et aLlbOOOD we find that 
the present burst would last at least 100 Myr at its current 
rate. These time scales show that the gas consumption in the 
galaxy is currently very rapid but its future development is 
highly uncertain. 

Finally, the similarity between this galaxy and the high- 
redshift LBGs provides an opportunity to constrain the star- 
forming mechanisms of the earliest galaxies, when the metal- 
licity was low and the environments were unrelaxed and 
chaotic. 

This work, together with the analysis done by 
lOstlin et al.l (|2003l ) and two future publications on the BCGs 
ESO 185-IG03 and Mrk 930, will be part of a statistical anal- 
ysis to constrain and define the major properties of the SSC 
populations in BCGs. 

The paper is organized as follows. In section 2 and 3 we 
describe in detail our data reduction and cluster detection 
method. Section 4 focuses on the estimation of the clusters' 
physical properties, in particular ages, masses and extinc- 
tion. In this section we also assess the evidence for fiux ex- 
cesses at red wavelengths. Section 5 we examine the possible 
causes of the red excess. In section 6 we report on the prop- 
erties of the host galaxy as showed by the cluster analysis. 
In section 7 we collect our main conclusions. 

The luminosity distance and distance modulus of Haro 
11 are estimated to be 82.3 Mpc and (m-M)=34.58 magj, 
respectively, assuming a cosmology of Ho = 73 Km s~^ 
Mpc"\ Qm = 0.27, and Qa = 0.73 throughout. 



2 DATA SAMPLE AND REDUCTIONS 

A summary of the full data set used in the present analy- 
sifl is given in Table [T] We have used archival HST/ACS 
imaging of Haro 11, with observations from the FUV to 
the optical (~ SOOOA) wavebands. We have supplemented 
these data with new[j HST/WFPC2 broadband imaging in 
V (F606W), I (F814W) and H (NIC 3/F160W; NICMOS 
near infrared NIR). To sample further into the infrared we 
have used the NaCo adaptive optics (AO) imager in the Ks 
banqj. Unless otherwise specified, we use the Vega magni- 
tude system throughout the paper. The zero points in Vega 
magnitudes are given in Table [l] for all the filters. 



task (|Koekemoer et all (|2002l 'l: iFruchter fc Hookl (|2002l ')') in 
PyRAF/STSDA^ and rescaled to the HRC pixel scale, that is 
0.025" /pi. The final science images are rotated to North-up. 
The short i/i?C F814W exposure was only used for estimat- 
ing the flux of the two brightest sources, in the center of 
knots B and C (Figure [T]). Both knot centers were saturated 
in the deep WFPC2 F606W and F814W exposures. 



2.2 WFPC2 data reduction 

Due to the small size (major axes ~ 13.0") of the star- 
burst regions into the galaxy, the target could easily fit into 
the high-resolution (pixel scale is 0.045") Planetary Camera 
(PC) which has a FOV of 35" x 35" . In order to refine the 
WFPC2 resolution we used sub-pixel dithering and com- 
bining patterns. The data reduction were performed with 
MULTIDRIZZLE. We took into account sky corrections dur- 
ing the photometric analysis rather than performing sky- 
subtraction with MULTIDRIZZLE. Before combining the final 
image, we resized the pixel scale to 0.025", applying a drop 
size of 0.8 for both the filters, and rotated the final image 
so that North was up. 



2.2.1 The WFPC2 point spread function 

Because we found no isolated bright objects either in the 
galaxy or in the PC field, we constructed a filter-dependent 
point spread function (PSF) using archival data. We re- 
trieved WFPC2/¥?>14m images of the NGC-104-F field of 
the globular cluster 47 Tuc (# GO 8160, PI I. King). This 
field was chosen because it is not crowded and has dither- 
ing path and exposure time comparable to our data. We 
built the F606W PSF in a similar fashione using data from 
another uncrowded stellar field (# GO 8090, PI S. Caser- 
tano). We reduced the images applying the same method, 
taking into account only the PC field of the mosaic. Iso- 
lated sources in the field were selected by eye, and encircled 
energy distributions for each of them were inspected in or- 
der to avoid sources with irregular PSF shapes. The tasks 
PSTSELECT and PSF were used to produce the final PSF. 
These two constructed PSFs were used to perform a com- 
pleteness test (see section 3.2). 



2.1 ACS data reduction 

A full descripti on of the d a ta re duction of the archival 
data is given in lOstlin et al.l ([20091). The final ACS images 
have all been drizzled and aligned using the MULTIDRIZZLE 



1 value taken from NASA/IPAC Extragalactic Database (NED) 
■^ Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science Insti- 
tute, which is operated by the Association of Universities for Re- 
search in Astronomy, Inc., under NASA contract NAS 5-26555. 
3 Associated with programs # GO 9470 (PI D. Kunth) and # 
GO 10575 (PI G. Ostlin). 

" Associated with programs # GO 10902, PI: G. Ostlin. 
* Based on observations made with ESO Telescopes at the 
Paranal Observatory under programme IDs 079.B-0585 and 
081.B-0234. The observations were carried out the July 24**", 
2008. 



2.3 NICMOS data reduction 

The NIC3 camera was chosen because of the size of the field 
of view, 51.2" x 51.2". To sub-sample and improve the pixel 
scale resolution (for NIC3, the plate scale is ~ 0.2" /pa;) a 
double NIC-SPIRAL-DITH path with 3 dithered pointings 
each was applied in order. The single science frames showed 
significant differences in the sky background levels and pres- 
ence of bias residuals, which we corrected before running 
MULTIDRIZZLE. In order to get the best combination of pixel 
sub-sampling and pixel drop size values we compared the 
flux growth curves of some bright point-like objects in the 
flnal science frame with a reference model PSF produced 



" STSDAS and PyRAF are products of the Space Telescope Sci- 
ence Institute, which is operated by AURA for NASA 
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Table 1. Haro 11: ff5T observations. Archival observations are associated with programs # GO 9470 (PI D. Kunth) and # GO 10575 
(PI G. Ostlin). The new observations were obtained in the program # GO 10902 (PI G. Ostlin). The zero point in Vega magnitude 
system, and the aperture correction, ac to the point source photometry are included for each filter. The last two columns show the 
number of objects detected with a S/N ^ 5, and the corresponding magnitude limit in each frame, (a) Shallower image used only to do 
photometry in the two central cluster regions saturated in the deeper exposure. 



Filter 


Instrument 


Exposure time 


ZP (mag) 


a.c(mag) 


N(cr s: 0.2) 


mag limit 


Archival obs. 


F140LP (FUV) 


ACS/SBC 


2700 s 


20.92 


-0.49±0.05 


113 


25.1 


F220W (NUV) 


ACS/HRC 


1513 s 


21.88 


-0.41±0.05 


80 


23.2 


F330W (U) 


ACS/HRC 


800 s 


22.91 


-0.35±0.05 


90 


23.7 


F435W (B) 


ACS/WFC 


680 s 


25.79 


-0.50±0.05 


178 


27.5 


F550M (V) 


ACS/WFC 


471 s 


24.88 


-0.45±0.05 


144 


26.2 


F814W (I) 


ACS/HRC 


100 s 


24.86 


-0.56±0.07 


(a) 


(a) 


New obs. 


F606W (R) 


WFPC2/PC 


4000 s 


22.89 


-0.61±0.09 


211 


27.5 


F814W (I) 


WFPC2/PC 


5000 s 


21.64 


-0.73±0.04 


211 


26.2 


F160W (H) 


NIC 3 


5000 s 


21.88 


-2.45±0.37 


71 


25.6 


K, 


VLT/NaCo 


1222 s 


23.32 


-2.15±0.64 


66 


20.5 




Figure 1. The starburst regions in Haro 11. The image is the result of the combination of the WFPC2 F814W filter in red, and the 
two ACS filters: F435W i n green and F 2 20W in blue. The orientation and the 1" length are indicated. The names of the three starburst 
regions are the same as in lKunth et al.l 1120031 ). 
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by TinyTim softwarqj. The resulting best set of input val- 
ues for the final combined image was found to be pixscale 
= 0.067" /px and pixfrac — 0.7. The pixel scale applied to 
the final frame is the diffraction limited value for the NIC3 
camera. The final science frame was rotated to North-up 
and calibrated in units of DN/s. We performed a number 
of tests of the calibration of the final IR frame: these are 
detailed in Appendix A. 

2.4 NaCo data reduction 

Ks imaging was carried out at the VLT UT4 telescope us- 
ing Naos-Conica (NaCo) AO instrument in natural guide 
star mode. The CONICA imager and spectrograph equipped 
with the AladdinS array and S27 camera (FOV of 28" x 28") 
was chosen. The chopping mode and dithering path were set 
and a total of 24 target images and 24 sky fields were taken, 
each with PIT = 9.4s averag ed on 13 NDITs. We used 
PyPELINE (JMicheva et al.ll2010l ) for bad-pixel masking, pair 
subtraction, flat field correction, source masking, sky back- 
ground subtraction and alignment. Before the final stack- 
ing, the quality of every frame was checked and those with 
blurred PSFs were excluded, in order to achieve the best pos- 
sible resolution. Frames with PSF full width at half maxi- 
mum (FWHM) > 4p2:(~ 0.108") were rejected. The final 
image was built from the median of the best 10 selected 
frames (corresponding to a total exposure time of 1222 s). 
A standard star was observed the same night, with the same 
instrumental set, and at comparable airmass values that for 
our target (the difference is < 0.08) so that the ZP value 
is affected by the same atmospheric extinction value as the 
galaxy. The pixel scale in the final frame was 0.027" /px. 



3 CLUSTER SELECTION AND 
PHOTOMETRY 

3.1 Cluster detection 



SExtractor (jBertin fc ArnoutsI Il996l ) was used to find 
sources in the deep R and I WFPC2 exposures. The star- 
burst in the galaxy is intense and many point-like sources lie 
in crowded regions and superimposed on the bright galac- 
tic background. Farther out of the main body, point-like 
sources can be still observed but now in empty regions with 
lower background values. In order to extract candidate clus- 
ters in all these different regions of the galaxy field we ran 
SExtractor twice on each frame, each time with different 
input parameters. To find sources in the outskirts of the 
galaxy, we used a detection threshold of 2a above the back- 
ground RMS noise, a minimum of 4 adjacent pixels and a 
minimum contrast for de-blending set to the default value, 
0.001. To detect sources into the starburst region we in- 
creased the detection threshold to 2.5cr, applied one of the 
mexhat filters available in SExtractor, and raised the value 
of the de-blending minimum contrast to 0.005. In this way 
we avoided large numbers of spurious detections while at the 
same time being able to separate blended sources in crowded 
regions. The two lists of coordinates were then merged into 
a unique catalogue, cleaned from multiple detections of the 

"^ http://www.stsci.edu/software/tinytim/ 
© 2010 RAS, MNRAS OOO.ITH231 




Figure 2. Colour-magnitude diagram for the 563 cluster candi- 
dates detected in both R and I filters. The black filled dots are 
the 211 sources with a S/N ^ 5 in both filters (see Table [TJ. 
Model evolutionary tra cks with metallicity Z = 0.004 are from 
IZackrisson et al.l 1 120011 ). The masses corresponding to each track 
are indicated. The thick squares show the model starting point 
at 1 Myr, the triangles the end at 14 Gyr. The arrow shows the 
direction of a correction for extinction of E{B — V) = 0.3. 



same object and from sources separated by less than the 
diffraction limit of the telescope. We fixed a center for the 
galaxy (a = 00 : 36 : 52 : 62.7, 5 = -33 : 33 : 17.49) and 
delimited detections to a region of radius of 13.75" around 
the center. Finally, detections in the two cleaned R and I 
band catalogues were cross-checked. All the sources which 
were not detected in both R and I catalogues were discarded. 
The final catalogue comprised 563 objects. 

For the two WFPC2 images, we performed a correction 
for CTE (charge transfer ef ficiency) effects on the observed 
fluxes using the algorithm of lDolphinI ([20091). Because of the 
complexity of the region, i.e. point like sources located in the 
extended luminous background of the galaxy, we chose not to 
use the median background value, as prescribed by Dolphin's 
formula. Instead we estimated a " local" background for each 
source, given by the median value of the pixel columns above 
the source, where few a percent of the total charge has prob- 
ably been trapped. 

We show the 563 cluster candidates in the colour- 
magnitude diagram (CMD) in Figure [2] 



3.2 Completeness limits 

To test the method we used to constrain the source catalogue 
and the detection limits, we performed a completeness anal- 
ysis for both R and I WFPC2 frames. A mock catalogue was 
made with objects positioned along a grid of 30 pixels per 
side (big enough not to increase crowding) . All the positions 
were located inside the circle defined previusly. Magnitudes 
were randomly assigned to each position, from 19.0 mag to 
29.0 mag with an interval of 0.5. For each filter, we created 
a synthetic frame with the sources using mk synth, a tas k 
of the in the public software toolkit Baolab (|LarserJll999l ). 
We supplied the PSF shape. The object detection was done 
in the same way as previously explained and no conditions 
were placed on the photometry. We then counted how many 
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Figure 3. Recovered fractions of sources as function of their 
magnitude for R and I images. The two continuous thick lines 
are the completeness limits for detection in filters F814W (black) 
and F606W (green) . The dotted (dashed) line show the recovered 
fractions inside (outside) the crowded starburst region, as defined 
in section [3. II The two red dotted horizontal lines show the 90% 
and 50% completeness levels. 



of the input objects were recovered, allowing a displacement 
not larger than 1.5 pixels. We repeated the procedure three 
times, shifting the grid by 10 x 10 pixels and re-assigning 
magnitude randomly to prevent introduction of systematics. 
The mean recovered fractions of objects as function of mag- 
nitude is shown in Figure [S] The completeness goes down to 
90 % at mag 27 in both filters, and only 50% of objects are 
recovered at ~ 28.0 mag. 

This is an averaged result: in reality, incompleteness 
due to crowding in the central active star-forming region is 
larger than in the outskirts of the galaxy. To assess by how 
much, we did a second completeness estimation, this time 
distinguishing between the two regions, with boundary ra- 
dius fixed at 6.5" from the center. The fraction of recovered 
objects inside the crowded region (dashed and dotted lines 
in Figure ^ begins to diminish in both R and I frames at 
between 25.0 mag and 25.5 mag, reaching 90% already at 
~ 26.0 mag and dropping to 70% at 27.0 mag. On the other 
hand, the recovered fraction of objects outside the crowded 
region gives 100% completeness for detections as faint as 
27.0 mag. As expected, blending affects the detections at 
the faint end of the cluster luminosity function in the active 
starbursting region. 



3.3 Photometry and final catalogue 

Aperture photometry was performed with the task PHOT, in 
all the available filters, using the catalogue previously deter- 
mined. The Iraf tasks GEDMAP and GEDXYTRAN were used to 
transform the coordinates from the WFPC2 system to the 
corresponding ACS, NICMOS, and NaCo/S27 catalogues. 
Growth curves of several objects in different filters were in- 
vestigated to decide the size of the aperture radius. Finally, 
we carried out photometry in all frames using a fixed aper- 
ture radius of a^ = 0.1" and estimating the local sky back- 
ground in an annulus of internal radius Sr = 0.125" and 
width 0.05". In order to to avoid photometric uncertainties 



affecting the age and extinction estimates, we rejected for 
each filter data points with magnitude errors Om > 0.2. In 
Table [T] we show the corresponding magnitude limits and 
the number of objects in each frame after this selection. For 
the final catalogue, we retained only objects detected in at 
least 3 filters. In total it contains 198 objectqj. Galactic ex- 
tinction corrections in all the filters was applied. Only 17% 
(34/198) of sources had detection in 3 fihers, 21% in 4, 12% 
in 5, 11% in 6, 13% in 7, 12% in 8, and 14% (28/198) in 
the 9 available filters. In total, 83% of the objects in the 
catalogue had detections in at least 4 filters. 

Aperture corrections and associated errors were esti- 
mated for each filter frame and added to the fluxes and 
photometric errors which we previously determined. They 
are given in Table [T] We used TinyTim PSF simulator to es- 
timate the aperture corrections values for ACS and WFPC2 
frames. A different method was applied to the two IR filter 
frames, H and Kg. The H frame has a wider PSF and lower 
pixel scale resolution. For these reasons, we decided to es- 
timate the aperture correction directly from isolated bright 
star-like objects in the frame. In the large FOV of the NIC3 
camera we could identify some point like objects. Moreover 
an identical set of NICMOS data was available for other two 
BCGs, Mrk 930 and ESO 185- 1 3 (par t of the same project 
# GO10902, see lAdamo et all (|2009l ): resuhs wiU be pub- 
lished in a following up paper, Adamo et al., in prep.) and 
also in those frames, reduced in the same way, few star-like 
objects were available. The aperture used to do photometry 
(0.1") was quite small and corresponded to half a size of the 
original NIC3 pixel scale (0.2"), making the error due to the 
aperture correction large. On the other hand, using bigger 
radius apertures would introduce blending of sources. 

The aperture correction for the Ks image was derived 
from 3 bright isolated point-like objects in the frame. In this 
case the problem was not the poor pixel resolution scale, but 
the effect of the AO technique on the PSF shape. The flux 
of the source is split between the diffraction-limited gaus- 
sian PSF core and the extended wings, making the aperture 
correction PSF-dependent. As expected, the aperture cor- 
rection value is high ac — —2.15 ± 0.64, because a consider- 
able part of the flux is contained in the PSF wings and the 
associated error is large. 



4 SED FITTING AND CLUSTER ANALYSIS 

4.1 Models 

For our SED fit s of th e target clusters, we have used the 
IZackrisson et al.l (|200ll . hereafter ZOl) spectral synthesis 
model. This model predicts the combined SED of both stars 
and photoionized gas. Often neglected, the gas component 
can actually have a pronounced impact on t he broadband 



fluxes of very young stellar popula t ions ( e.g. Kriiger et al 



2008 



19951: lAnders fc Fritze-AlvenslebenI l2003l : IZackrisson et al 



The ZOl model, which includes pre-main sequence evo- 
lution and a stochastic treatment of horizontal branch 



The complete photometric catalogue is available at 
|http://www. astro. su.sc /~adamo or on request from the au- 
thors. 
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morphologies at low metalli cities, is bas e d on syn- 
thetic stellar atmospheres by iLeieune et alj (| 19981 ) and 
IClegg fc Middlemasd 11987]), together with stellar evolution- 
ary tracks mainly from the Geneva group. The gas contin- 
uum and emission lines are predic ted using the photo ioni- 
sation code Cloudy, version 90.05 (|Ferland et al.|[l99a ). At 
each time step, the spectral energy distribution of the in- 
tegrated stellar population is fed into Cloudy, which then 
provides a realistic evolution of the emission-line ratios over 
time. The spectra predicted by the models are redshifted to 
match the target galaxies and convolved with the transmis- 
sion profiles of the filters used. 

In the ZOl model, we assume a stellar metallicity of 
^stars = 0.004, a Salpeter initial mass function {dN/dM oc 
j^-2.35^ throughout the mass range 0.08-120 Mq and an in- 
stantaneous burst (i.e. single-age) population. A conversion 
factor is appli ed to correct the estimated cluster masses to a 
Kroupa IMF (iKroupalbOOll ) in the mass range 0.01-120 Mq. 
The properties of the nebular component are determined by 
the adopted gaseous metallicity Z^as, the hydrogen number 
density n(H), the gas filling factor / and the gas mass Mgas 
available for star formation. We adopt n(H) = 100 cm~^ and 
/ = 0.01 (values typical for Hll regions), a gas metallicity 



identical to the stellar one {Zg, 



and a gas mass 



of Mgas = 10 M(7). Finally, we hav e used the spectral syn- 
thesis model of lMarigo et all (|2008l . hereafter M08) to test 
the robustness of some of our conclusions. This model fea- 
tures a more sophisticated treatment of thermally pulsating 
asymptotic giant branch (AGB) stars (which mainly affects 
the interpretation of near-infrared data), but does not in- 
clude nebular emission. For this reason it i s not suitable for 
the a nalysis of clusters at ages ^10 Myr (jZackrisson et al.l 
[200I). 



4.2 Constraining the physical properties of the 

SCs 

4-2.1 Chi square method 

We used a least-squares fit to estimate for each cluster the 
best fit to ou r models (similar to other previous works, see 
lAnders et al.ll2004-, : Bik ct al.; 20031 ). Internal extinction was 
treated as a free parameter, and we used the Calzetti ex- 
tinction law (|Calzetti et al.ll2000l ). allowing extinction vary- 
ing from E{B - V) = 0.0 to 3.0 with a step of 0.01. The 
Q-valuJQ and the reduced Xr = X^ /'^ were calculated. The 
u = A'' — m, is the number the degrees of freedom. In our case, 
for each cluster, the number of detections in different filters, 
A ^ 3, while the number of parameters m = 3 (the age, 
mass and extinction of the cluster). The Q-val ue is a sta- 
tistic al measurement of the goodness of the fit (|Press et al.l 
[1993). Q and Xr were estimated only for clusters with at 
least detection in 4 different filters (A > 3). 

Among our initial model fits, we noticed several cases 
for which Xr >> 1 f^nd Q < 10~^. For those cases, the in- 
spection of the residuals (displacements between the mod- 



^ also known as the x^ probability function, see iPress et al.l 
1119921 ). Best-fit models with values of Q which exceed 0.1, are 
good representative of the observed data; values of Q > 0.001 
can be considered still acceptable; while models with lower values 
are unlikely. 



eled and the observed integrated photometry of the cluster. 
Am = TTiinod — rnobs), revealed an unusual positive offset, or 
"excess", of the observed fluxes in H, Ks, and I bands. 

We performed two further sets of x^ fits, this time ex- 
cluding the NIR data. The first of these we refer to as UV- 
UBVR, where only the 6 bluer filters were included in the 
fit. The second set included 7 filters, FUV, NUV, U, B, V, 
R, and I, indicated as UV-UBVRI. The initial fit, includ- 
ing all the available data from UV to near-IR, we denote 
UV-UBVRIHKs- We maintained the requirement of at least 
detections in 3 filters for each SED fit. This restriction meant 
that we could fit 190 objects when the H and Ks bands were 
excluded in the UV-UBVRI fit. Of these, 35% have detec- 
tions in all the 7 included bands, while 22% have detections 
in only 3 filters. When the I band was excluded, the num- 
ber of objects fitted decreased to 149 (23% with at least 3 
detections and 45% with detections in all the allowed filters) . 

4-2.2 Investigating the excess at near-IR wavelengths 

In this section we present "observational" evidence for the 
red excess and at the same time probe that the fit to the 
UV-blue optical integrated fluxes of the clusters produces 
robust constraints of the cluster properties. 

In Figure |3] we show an example of the three differ- 
ent fits to the integrated SED data points for four different 
clusters with different physical properties (a table with inte- 
grated photometry of these clusters is presented in Appendix 
B). The quantities plotted have been converted from Vega 
magnitudes to physical fiux units, erg s~^cm~'^A~^. Inte- 
grated photometric data points with associated error bars 
are plotted together with the best synthetic data points ob- 
tained by the three different types of fit: UV-UBVRIHKs 
(black squares), [/F-f/B FT?/ (red triangles) and UV-UBVR 
(blue diamonds). The underlying spectrum represent the 
best model. The broad F606W filter transmits the promi- 
nent Ha emission line. So in cases where we see a pronounced 
jump between fiuxes in the narrow F550M and the broad 
F606W filter we expect young (^8.5 Myr) stellar complexes 
(e.g. see clusters #43 and #113 in Figure |3|. The displace- 
ment between these two data points is negligible already 
around 10 Myr (e.g. cluster #6 in Figure Q. Finally, the 
data points sampling the U-UV part of the spectrum are 
extremely sensitive to the internal extinction of the cluster. 
In Table [5] we list the statistical results and best fit param- 
eters obtained with all three sets of fits. 

The fit to cluster # 113 is quite good all the way from 
the UV to the IR (inside the error bars). Looking at the 
Table (2] we see that removing filters from the fit produces 
no noticeable effect. The jump between V and R clearly 
confirms the young age of this object. The non-detection 
in the UV part of the spectrum is in agreement with the 
high extinction value. The properties of 83 clusters with 
similar SED fits to this cluster were estimated by the UV- 
UBVRIHKs fit. 

The UV-UBVRIHKs fit to the cluster #87 clearly failed 
to fit the integrated IR flux in Ks , sitting above the best-fit 
model with an excess of AraKs ~ 2.34 mag. The large error 
associated with the observed flux in Ks produces a smaller 
effect on the Xr and Q values of the UV-UBVRIHKs fit. 
Taking into account this factor, even small differences in 
Xr and Q, as observed here between the UV-UBVRIHKs 
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Figure 4. SED analysis of 4 clusters. The filled black points indicated the observed photometric values of each cluster. The integrated 
model fluxes are labelled with different symbols for each set of fit as indicated in the plots. We plotted the spectrum of the final best 
fitting model with a colour which corresponds to the type of fit. 



and UV-UBVRI fits, become important. In tliis case, we 
discarded the UV-UBVR fit solution because, as indicated 
in Table O the exclusion of tiie I band from the fit does 
not give any improvement. We found 26 clusters witii an 
excess at the IR wavelengtlis (H and Ka). For these objects, 
the final ages, masses and extinctions were derived from a 
UV-UBVRI fit. 

The fits to cluster #43 give even more dramatic re- 
sults. Neither the UV-UBVRIHKs nor the UV-UBVRI His 
are good (high Xr and high Q). The estimated age, 45 Myr, 
is impossible to reconcile with the jump observed between 
the narrow V and broad R filters, and the estimated mass 
is high. Excluding the I band photometry gave a better fit 
(see statistical values in the Table 2). The UV-UBVR fit to 
this cluster produces acceptable Xr a-nd Q, indicating a clus- 
ter only 3.5 Myr old, with mass ~ IO'^Mq. E(B-V) is 0.34, 
which corresponds to a visual extinction of Av = 1.24. The 
observed excesses compared to the best UV-UBVR model 
are 0.91, 2.01 and 2.97 mag in I, H and Ks, respectively. 
A similar case is cluster #6. Here, too, the fiux excess is 
not confined to the NIR part of the spectrum, but also ex- 
tends into the I-band, if less dramatically than #43. In #6 
the V-R jump is also less pronounced, in agreement with its 
best estimated age, 8.5 Myr. We found a total of 78 objects 
affected by flux excess at wavelengths > 8000 A like clus- 



ters #6 and #43. For these clusters, we have determined 
physical properties from the UV-UBVR fit. 

After this first look, it is evident that while our SSP 
models can successfully model the optical and UV portions 
of the observed cluster SEDs, they are unable to reproduce 
the NIR photometry of many of the clusters. This suggests 
that the models we have used are missing a source of emis- 
sion which is bright in the near-infrared. 

In Figure[5]we show how the three different fit types af- 
fect the residuals in the SED fits in all filters, for targets with 
flux excess in the H band. The black hatched histograms are 
the residuals from the UV-UBVRIHKs fit relative to the ob- 
served data points. Am = TTimod — Jriobs, for all the clusters. 
The residuals in H and Ka are systematically offset toward 
positive Am. In all the other filters the spreads are large, but 
the residuals are evenly distributed around zero. When the 
UV-UBVRI And UV-UBVR fits were performed, the spread 
in the residuals of the excluded filters (red thick histograms 
for F160W, Ks, and F814W ) were even larger. On the other 
hand, the improvements in the distributions for R and blue- 
ward filters are clear. The red histogram of the F814W filter 
has a secondary peak around zero. This peak is produced 
by the clusters which do not have any flux excess in I band 
and for this reason fltted by the UV-UBVRI fit. 

In Figure [G] we show four different colour-colour dia- 
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Figure 6. Colour-colour diagrams for the cluster samples. The filled black dots are clusters without a red excess, fitted on UV- 
UBVRIHKs- The open red triangles are clusters with flux excess in H and K^ and fitted on UV-UBVRI, and the blue diamonds denote 
clusters with excess in all of I, H and Kg, fitted only on UV-UBVR. The thick black line shows the predicted stellar evolutionary model 
according to ZOl. The black dashed line is the corresponding M08 model, with Z = 0.004. The black triangles indicate the youngest ages 
included in each model, 0.55 Myr for ZOl and 4 Myr for M08. The red crosses mark the position in the tracks corresponding to 10 Myr 
and the black squares the end of the tracks at 14 Gyr. In each plot the mean photometric error associated with each object is shown. 
The arrows indicate direction and displacement produced by an extinction correction of E[B — V) = 0.3. The number of objects varies 
from plot to plot due to different detection limits in each filter. 



grams for the cluster sample together with the evolutionary 
tracks by ZOl for diflferent filter combinations. M08 evolu- 
tionary tracks with the same metallicity {Z — 0.004) are 
also shown. As pointed out in the previous section, the in- 
clusion of nebular continuum and emission treatment re- 
ally do make a fundamental difference in the models at the 
youngest ages, and at all the wavelengths (see differences 
between the ZOl tracks and those of M08). In particular the 
broad filters (F435W, F606W, F814W, F160W) are more af- 
fected, with colours changing drastically in the first 10 Myr 
(see F606W-F814W, F435W-F606W, F160W-F606W). 
Finally, the colour diagrams which include IR filters also 
show a difference at older ages (~ 1 Gyr) between the two 
evolutionary models which can be attributed to the better 
treatment of the AGB-phase by M08. In the next section we 
will show, however, that because of the youth of the clusters 
in Haro 11, this difference actually only affects a very small 
fraction of the sample. 

The differences in the number of objects in the colour- 
colour diagrams are due to the different detections limits 
in each filter. We represent the three cluster subsamples by 
coloured symbols corresponding to the flux excess. Black 



flUed dots are clusters which do not show any flux excess, 
red open triangles clusters with observed flux excess in H 
and Ks; blue open diamonds are clusters with an additional 
observed displacement in I, and cluster with excess in I and 
none detections in the IR bands. In the optical and UV 
colours (the upper panels of Figure [6| , we see no difference 
in the distributions of the clusters with and without flux ex- 
cesses in I band. In these panels, many of the clusters can be 
explained as a young cluster population affected by differ- 
ent values of extinction. Few clusters, affected by red excess, 
show very red F606W-F814W colour {FQ06W - F814VF > 
1.0 mag) and are displaced relative to the tracks. These clus- 
ters and many of the blue open diamonds which show red 
F606W-F814W colours {F606W - F8UW ^ 0.5 mag) have 
UV colours (left column, upper row) consistent with stellar 
populations younger than 10 Myr (F140LP - F606W ~ -2 
mag). 

Including the IR wavebands (Figure [O] lower panels) 
in the colors, changes the situation drastically. Clusters af- 
fected by red excess are now displaced from the model tracks 
by up to 4 magnitudes. Such offsets are impossible to rec- 
oncile with the estimated values of extinction. Moreover, 



© 2010 RAS, MNRAS 000.[TH23l 



10 A. Adamo et al. 



Table 2. Final outputs given by the three sets of SED fits as 
shown in Figure |4] In bold we show the final age, mass and ex- 
tinction assigned to the clusters. 



id 



UV-UBVRIHKs 
Q Myr IO^Mq E{B-V) 



113 


0.7 


0.72 


5.5 


1.4 


0.44 


87 


1.65 


0.25 


10.5 


1.3 


0.0 


43 


13.6 


< 10-^ 


45.0 


16.3 


0.23 


6 


3.28 


0.006 


12.5 


0.8 


0.0 



UV-UBVRI 



id 



Myr IO^Mq E(B-V) 



113 


0.45 


0.5 


5.5 


1.4 


0.44 


87 


1.1 


0.34 


10.5 


1.3 


0.0 


43 


13.1 


< 10-^ 


45 


16.0 


0.23 
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2.6 


0.008 


12.5 


0.8 


0.0 
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Myr lO^M© E(B-V) 
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0.5 


87 


0.8 


0.31 
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1.2 


0.0 


43 


1.4 


0.07 
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0.41 
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0.01 








Figure 5. Residuals (Am = rrimod ~ "lobs) of the x^ SED fit- 
ting in all the available filters for sources with detected excess 
in H. The hatched black histograms show all these clusters; the 
histograms with red borders arc the residual distribution for the 
same targets but where I (if the flux excess is present), H and Ks 
have been excluded from the SED fltting. The vertical thick dot- 
ted lines indicates zero residuals, i.e. a match between model and 
observation. The number of objects differs between plots because 
of the different detection thresholds in each band. 



several of the sources with red F606W-F814W colour and 
strong UV emission, which we expect to be very young, have 
IR colours which lie close to the oldest ages of the model 
tracks at ~ 14 Gyr. This is again impossible to explain with 
the current models. 

We analyze possible mechanisms producing the excess 
at the redder wavebands in section [S] In the next section we 
will focus on the derived physical properties of the clusters. 




Figure 7. Final age distribution of the star cluster population 
(black hatched histogram). The dotted over-plotted distribution 
shows the age distribution if the SED is fitted over the whole UV- 
UBVRIHKs range, without taking the red excess into account. 



4-2.3 Final estimated ages, masses and extinctions 

In Figures [T] [HI and [5] we show our adopted age, mass, and 
E{B-'V) distributions for the total cluster sample, with 187 
objects (black hatched histograms). In this final sample, only 
in 21% of the cases we determined cluster properties using 
3 filters, while the 53% had detection in at least 5 or more 
bands. The cluster properties of the targets affected by a 
red excess have been obtained from the fit to the blue-side 
of the SED spectra. For comparison we over-plotted the cor- 
responding distributions (dotted histograms) produced by a 
UV-UBVRIHKs fit to the entire sample of clusters. The red 
excess clearly introduces a systematic offset in the estimates 
of all the derived parameters. When all the filters are in- 
cluded in the fit, clusters with a red excess appear to be older 
and more massive than their UV-optical photometry alone 
would indicate. What happens is that when we attempt to 
fit the NIR filters affected by flux excess, the best model 
choices are old stellar populations in which lower-mass stars 
dominate the light. This affects the mass estimation, because 
the M/L ratio increases with age: for the same luminosity, 
older clusters are also more massive. For example, cluster # 
43 (section 14.2.21 and Figure [Jjl shows that even the excess 
in I band (at ~ 8000 A) has important consequences for the 
estimation of the cluster ages and masses. 

So far, the fllter F814W is one of the most used in 
the photometric analysis of the cluster populations in extra- 
galactic environments. Because the flux excess in this filter 
is not as dramatic as in the IR filters, it may go unnoticed 
when IR data are not available. Indeed, this affe cted our own 
initia l analysis of Haro 11 's cluster population (|Haves et al.l 
[20071). In that paper, the cluster age distribution showed the 
same double peak as we have shown here where no correction 
for red excess is performed (dotted histogram in Figure [7}. 
Extinction estimates, on the other hand, tend to be more 
sensitive to the UV-U fluxes and are therefore not affected 
by red excess. 

The cluster population in Haro 11 appears to be very 
young (Figure [TJ . The present starburst phase started less 
than 40 Myr ago with a peak at the extremely young age 
of 3.5 Myr. Only 16 clusters have estimated ages older then 
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Figure 8. Final mass distribution (blacli hatched histogram) of 
the star cluster population. The dotted over-plotted distribution 
shows the masses obtained for the clusters if the SED is fitted 
over the whole UV-UBVRIHKs range. 




Figure 9. Final extinction distribution (black hatched his- 
togram) of the star cluster population. The dotted over-plotted 
distribution shows the extinctions obtained for the clusters if the 
SED is fitted over the whole UV-UBVRIHKs range. 



40 Myr and of these, half have an age > 100 Myr. The age 
distribution is agreement with the esti mated burst duration , 
Tb = Mburat/SF-R = 35 Myr, found by lOstlin et all (|200ll ). 
where Mburst is the total mass contained in the burst and 
SFR. They estimated an upper limit age for the present 
burst of 35 Myr, taking the SFR from the total Ha lumi- 
nosity, uncorrected for extinction. The age peak at 3.5 Myr 
is consistent with the young stella r population (age < 5 
Myr) which iGalametz et al.l (|2009l ) identified as responsi- 
ble for the hard radiation field inferred from the observed 
ratio of [Ne in]/[Ne n]= 3.2 in the galaxy's integrated mid- 
IR Spitzer Infrared Spectrograph spectrum. With about 130 
clusters younger than 10 Myr, Haro 11 represents a unique 
opportunity to investigate the cluster formation process and 
evolution at t he youngest ph ase i n starburst s ystems. In re- 
cent reviews, iLarserJ (|2010l ) and iLadal (|201(J ) pointed out 
that clusters younger then ~ 5 Myr are still partially em- 
bedded and affected by several magnitudes of internal ex- 
tinction due to the remains of the giant molecular clouds 
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Figure 10. Extinction versus cluster age. The small blue asterisks 
are the extinction as function of the age for each cluster. The black 
open diamonds indicate the mean value at each age step in the 
model, as long as at least three objects are contained in the bin. 
The associated error bars are the la uncertainties. 



in which they formed. In Haro 11, the estimated extinction 
distribution is quite wide with a peak at E{B — V) ^ 0.0 
and an extended tail down to at least 0.6 mag (Figure [O] 
corresponding to 0.0 ^ Ay ^ 2.2 mag). The few cases with 
E{B -V)^ 0.8 mag have 2.9 < Av < 3.7 mag. 

In Figure \TU\ we show how the derived extinction 
changes as function of cluster age. For ages between 1 and 
3.5 Myr the extinction range is wide. Almost all the clus- 
ters with E{B — V) > 0.4 (Av > 1.5) mag are located in 
this age range, but the number of clusters with lower ex- 
tinction is not negligible. This observed spread shows the 
complexity of the cluster formation process and the grad- 
ual dissolution of the parent clouds. A correlation might be 
expected between the extinction of clusters younger than 5 
Myr and their masses, but we found no clear trend. This lack 
of correlation could indicate that we have underestimated of 
the mass of partially embedded clusters. If this is the case, 
then we should consider the estimated masses at very young 
ages as a lower limit to the real values. At older cluster ages 
(> 5Myr), we can see in Figure [TOl that extinction tends 
to decrease as function of the age, although the deviations 
are large. Moreover, previous studies have shown that the 
galactic environm ent is quite complex (jBergvall Ostlirj|2002l . 
iHaves et al.ll2007l ). and the mean extinction in the galaxy is 
patchy. Some regions are apparently free from extinction 
while in others thick dusty clouds and filaments are appar- 
ent (see Figure[l]). Finally, we noticed that similar trends be- 
tween cluster ages and ex tinctions have also b een observed 
in t he Antennae sy stem (jMengel et al.l |2005| ) and in M51 
( Bastian et al.l |2005| ) . 

The clusters are quite massive. More than half of the 
targets have masses between 10* and 10^ Mq. In Haro 11 
we have found around 60 estimated SSCs (clusters with the 
masses > 10^ Mq), seven of which have masses larger than 
W'^Mq. The two knots, B and C, have masses ~ IO'^Mq. 
The fit for knot C is poor, however, and its estimated age, 
mass and extinction have not been included in the Figures [3 
[8] and (9] We present a more detailed analysis of these two 
regions below. Only ~ 20 clusters have masses below 10* M©, 
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Figure 11. Ages versus masses for the star cluster population. 
Estimates from UV-UBVRIHKs fits are shown as filled black 
dots, while estimates from UV-UBVRI fits and UV-UBVR fits 
are shown as red triangles and blue diamonds, respectively. The 
lines are the expected detection limits in various filters (see the 
labels) according to our stellar evolutionary models. 



though this small number is affected by the detection limits 
we have imposed on the available data. 

In Figure [11] we show estimated masses as a function of 
our derived cluster age s. The lines indicate the predicted de- 
tection limits using the lZackrisson et al.l (|200ll ) stellar evolu- 
tionary model, for the corresponding magnitude limits (see 
Table 1) as indicated in the plot. Detections below these 
thresholds are highly unlikely. All the clusters in the plot 
have been detected in R and I band (used to make the ini- 
tial catalogue; section 13. 3p . and at least in one other fil- 
ter. While the clusters affected by red excess have a wide 
mass range, we noticed that all the most massive clusters 
(M ^ 5 X 10^ Mq) younger than 40 Myr, appear to be af- 
fected by at least some flux excess in the NIR. We detected 
no clusters older than 40 Myr with any red excess. Six ob- 
jects with masses around 10^ Af© and ages greater than 500 
Myr have properties similar to more evolved globular clus- 
ters, as in the Milky Way. 

Curiously, even though restricted detection limits allow 
us to detect very young clusters (1 — 3 Myr) with masses 
around 10'^ — 10* Mq, we see almost none in this range of 
age and mass. This peculiar phenomenon could be inter- 
preted in different ways. (1) Low mass clusters < 10'' Myr 
are still embedded in their parental clouds so they are not 
"visible" in optical wavebands. (2) At the distance of Haro 
11, sources may be blended. Indeed, we expect cluster s to 
form in complexes rather than singly (|Elmegreenll2010l l. If 
this is the case, some of the more massive observed clusters 
may in reality be clumps of several low-mass clusters. (3) 
The model, in particular its handling of the nebular com- 
ponent, has limitations which may lead to objects at 1 — 3 
Myr having their ages overestimated. We will return to these 
three points in section 16.11 For the moment, we note that 
none of the three explanations excludes the others and all 
three could account for the apparent absence of very young 
low mass clusters. 

Finally, we have investigated the two most active star 
forming regions in the galaxy, knots B and C (Figure [T]). 
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Figure 12. SED analysis of the two knots, C (#1, above) and B 
(#2). 



We performed aperture photometry on the two central re- 
gions in the same way done for the point-like clusters even 
though these two regions appear extended. Knot C is more 
reminiscent of a nuclear cluster than a point-like SC. Dusty 
fllaments ap pear at the t wo sid es of the central region in a 
cone shape. iHaves et al.l (|2007l ) pointed out that knot C is 
the brightest source in the UV, with strong Lya emission 
from its centre, possible indicating that star formation in 
knot C is still ongoing. Knot B appears obscured by fila- 
mentary dust clouds passing through its central region. It 
is not clear from our images whether this central region is 
composed of a single massive object or several clusters in a 
clump similar to the ones obse r ved in the nearby starburs t 
galaxy M 82 (|Melo et all l2005l : IWestmoguette et all 120071 ) . 
Knot B is the most massive cluster-like object after C. Both 
knots appear affected by the red excess (see Figure I12|l and 
are very massive objects. The best fits were obtained with a 
UV-UBVRI at for both clusters. The physical properties of 
B and C are summarized in Table [121 As showed by the Xr 
value, the fit to knot C is poor so its properties have to be 
considered as an approximation. We have obtained spectra 
of these knots with VLT/X-shooter which will be the sub- 
ject of a future paper (Gumming et al., in prep.). Here we 
exclude knot C from our final analysis of the cluster popu- 
lation properties. 
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Table 3. Ages, masses, and extinctions for the two star-forming 
regions B and C. 



id 



age(Myr) mass(M0) E{B - V) 



B 

C 



3.3 

36.0 



3.5 
9.5 



8.35 X 10*' 0.38 

1.36 X 10^ 0.06 
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Figure 13. Recovered age distributions using Calzetti (green 
hatched histogram) and LMC (dotted histogram) extinction laws. 
In the inset, we compare for each cluster the age estimated with 
the two extinction laws. The different symbols are the same as in 
Figure 11. 



4.3 Can a different extinction law give different 
results? 

Could a different extinction law have an impact on our esti- 
mates of the physical properties of the clusters? To answer 
this question we performed a second x^ fitting run, this time 
assu ming a L arge Mag ellanic Cloud (LMC) dust extinction 
law (iFitzpatrick 1986: iMisselt etal.l [19981 . Both the LMC 
and Calzetti ( Calzetti et al.ll2000l ) laws have been obtained 
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Figure 14. Recovered mass distributions if a Calzetti (green 
hatched histogram) or an LMC (dotted histogram) extinction law 
is used. The inset shows a direct comparison is shown between the 
recovered masses for the two different fitting runs. The different 
symbols are the same as in Figure 11. 
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Figure 15. Recovered extinction distributions using Calzetti 
(green hatched histogram) and LMC (dotted histogram) extinc- 
tion laws. The inset shows for each cluster the corresponding 
visual extinction. Ay recovered in the two cases. The different 
symbols are the same as in Figure 11. 



from different environme nts and physical con dictions. The 
Calzetti extinction law (jCalzetti et al.l 120001 ) was derived 
from the integrated spectra of several starhurst galaxies, so 
it describes the total extinction affecting a galactic environ- 
ment and produced by different mechanisms. Since, dust is in 
general mixed in galaxies, different optical depths are probed 
as a function of wavelength, and the emerging 'Calzetti law' 
results from the combination of this effect and the under- 
lying extinction law. On the other hand, the LMC law was 
determined from the spectra of stars located inside and near 
the 30 Dor cluster and reproduces the local extinction pro- 
duced by a molecular dust cloud on single stars, in which 
case the dust can be accurately treated as a screen. The 
situtation in Haroll may, at this distance, be intermediate 
to these two cases, and most results presented assume the 
Calzetti law. However it was important to assess whether the 
LMC law produces different results. Comparing the outputs 
of the two runs, we noticed the following. (1) The values of 
Xr obtained with the LMC law are similar to the Calzetti 
law results and show no improvement in the fit. (2) The clus- 
ters still display pronounced red excesses. In Figures 1151 1141 
and[T3]we have plotted the recovered extinction, mass, and 
age distributions obtained using the Calzetti (green hatched 
histograms) and LMC (dotted histograms) laws. The use 
of a different extinction law does not have any significant 
impact on the age distribution of the clusters. Only a few 
clusters have discordant ages estimates, as we show in the 
inset in Figure [T51 On the other hand, about one third of the 
clusters fitted with an LMC extinction law appear to have 
higher estimated masses (Figure [14)). This trend in the es- 
timated mass can be explained by the recovered extinctions 
when the LMC extinction law is used. Comparing LMC and 
Calzetti E(B — V) distributions (Figure [15}, we see that a 
considerable number of objects with Calzetti E(B — V) in 
the range 0.1 — 0.6 appear distributed in a long tail toward 
higher extinction values. In particular, the quantitative dif- 
ference in visual extinction Ay can differ almost a factor of 
two, depending of which extinction law is used (see the inset 
in Figure [15)) . This in turn means that objects affected by 
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higher extinctions have higher estimated masses. While we 
cannot draw any definitive conclusion on the choice of the 
extinction law, it is clear that a considerable fraction of our 
recovered extinctions and masses can depend on the adopted 
extinction law. Such uncertainties should be taken into ac- 
count when cluster properties are estimated. Age estimates, 
on the other hand, appear to be robust. 



4.4 Testing the SED fit results with Monte Carlo 
simulations 

The outputs recovered by the SED fitting technique are con- 
nected to the stellar evolutionary models we used. In order 
to estimate the uncertainties in our recovered ages, masses 
and extinctions, we used Monte Carlo methods to simulate 
a fully sampled population of clusters which resembles the 
population we observe in Haro 11. The initial sample con- 
tains clusters with three different mass values (10'', 10^ and 
IQ^Mq) for each of 51 age step of the evolutionary model, 
and 3 extinctions {E{B — V) = 0.05, 0.3, 0.5), yielding a total 
of 3 X 51 X 3 = 459 clusters. 

For each cluster, we use the spectral synthesis models 
to calculate its apparent magnitude in all the 9 filters for 
the mass, age and E{B — V) in question. To these ideal 
magnitudes that perfectly match the model prediction, we 
add random photometric errors in accordance with the ob- 
served magnitude-error relation in each filter and assuming 
a gaussian error distribution. In this way we simulated 1000 
realizations of each cluster, the fitted them with our SED 
tool and compared the input to the output for 1000 hypo- 
thetical clusters. 

All the 1000 realizations of each cluster were fitted with 
the same Xr fitting program using an UV-UBVRIHKs fit, 
an UV-UBVRI fit, and a UV-UBVR fit in turn. The me- 
dian and quartiles of the distributions are considered to be 
the estimated output age and associated uncertainties, re- 
spectively, for each cluster. In Figure [TS] we show how the 
recovered age of an object and its uncertainty changes de- 
pending on the mass and extinction of the object, and the 
type of fit. From the initial cluster sample, combinations of 
age, mass and extinction which resulted in a cluster with 
apparent magnitude rriR > 27.5 mag was excluded from the 
analysis. Looking at the two panels of Figure [TS] we can see 
that this cut in luminosity appears as an incomplete sam- 
pling at older ages. In general, for the same mass, clusters 
with lower extinction are better fitted. The higher the clus- 
ter mass, the smaller the uncertainties at all ages and with 
respect to the reddening. For higher values of the reddening, 
input ages of ~ 10 Myr tend to lead to older output ages 
(and between 100 Myr and 1 Gyr) with a tendency to under- 
estimate the output ages. This last effect is less important 
for higher masses and for lower extinction. 

To test for systematic effects of our three different SED 
fits, we checked and compared the three sets of recovered 
ages for the same mass and extinction distributions. In Fig- 
ure [I6l we show two of the three SED fits. As expected, if 
fewer filters are available, the properties derived from the 
observed SED become more uncertain, and we see the quar- 
tile error bars get wider. The median of the distributions, on 
the other hand, are not affected except at two critical ages, 
around 10 Myr and 100 Myr. These age ranges correspond to 
particular features of the synthetic stellar population model. 
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Figure 17. Recovered age distributions (black hatched his- 
tograms) of the three sets of fits. Prom left to right, the re- 
covered age distributions for UV-UBVRIHKs, UV-UBVRI und 
UV-UBVR fits. The red histograms show the input ages of the 
simulated sample of clusters. 



There is a loop in the modeled colors at ages between 10 and 
40 Myr which makes it difficult to distinguish between dif- 
ferent ages steps. At intermediate ages, around 100 Myr, the 
tracks are parallel to the reddening vector (see Figure (6] up- 
per left panel) and clusters could be fitted equally well with 
younger ages and higher extinctions. We found no other ev- 
idence of systematics introduced by our analysis. This can 
be seen also in the plots in Figure [TT] in which we compare 
the input age distribution (histogram with thick red bor- 
der) of the simulated cluster population with the recovered 
one (black hatched histogram) derived from the three fit 
ranges. Once more we exclude clusters with m_R > 27.5 mag 
from the analysis. When all the filters are included in the 
UV-UBVRIHKs fit, the input and output distributions dif- 
fer little. With a UV-UBVRI fit the disagreement between 
the input and output distributions is more evident with a 
tendency to overpopulate the age bin around 10 Myr and 
60-80 Myr. The exclusion of the I band in the UV-UBVR fit 
moves the peak from 10 to 30 Myr, and the spurious peak 
a 60-80 Myr remains. In Figure [7] we see that the introduc- 
tion of a subsample of clusters analyzed with UV-UBVRI ot 
UV-UBVR fits produces a peak in the age distribution at 
3.5 Myr (log(age) = 6.3), populated by objects from older 
age bins, like 35 Myr. In the distributions of the simulated 
cluster sample we used the same age bin size as in Figure [T] 
None of the fits performed introduce any displacement to 
or spurious peak at 3.5 Myr. For this reason, we are quite 
confident about the robustness of our estimates. 



5 WHAT IS THE CAUSE OF THE RED 

EXCESS? 

Following our analysis of the cluster photometric and phys- 
ical properties, we address the question of what causes the 
excess of the observed fiuxes longward of 8000 Awith re- 
spect to the synthetic evolutionary models. Clusters with 
red excesses have also been observed in o ther galaxies 
IWhitmore fc Zhand[20o3 : ICresci et al.ll2005l : Reines et al. 
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Figure 16. Recovered ages as a function of the input ages for different values of masses and extinctions, as specified in the labels. Nine 
combinations of mass and extinction are shown for two of the fit ranges UV-UBVRIHKs (top), and UV-UBVR (bottom). The dashed 
lines mark where input and output ages are equal. The output ages are the median values of the 1000 realizations of each cluster and 
the error bars are the quartiles of the distribution. Simulated clusters with apparent magnitude mn > 27.5 mag have not been included 
in the analysis. This cut in magnitude is reflected in the plots as incompleteness of the full age sampling. 



2008a, b: iFernandez-Ontiveros et al.ll2009l l. They are typi- 
cally very young and bright at IR wavelengths, and appear to 
be subject to by significant dust extinction at optical wave- 
lengths. Their location in colour-colour diagrams tend to be 
offset toward both red V — I and V — K colors. An obser ved 
flux excess in I band was found bv lReines et al.l (|2008al ) in 
the young clusters of the galaxy NGC 4449. A subsequent 
spect roscopic analysis of two of these clusters (JReines et al.l 
120101 ) revealed that the flux excess could be explained if 



a contribution from nebular continuum and emission lines 
were taken into account. This explanation can be excluded in 
the case of the I-band excess in Haro ll's clusters, since our 
models already include both nebular continuum and emis- 
sion lines. 

Thanks to the large number of clusters with red excesses 
in Haro 11, we were able to check possible correlations be- 
tween the magnitude of the excess at A > SOOOA and the 
cluster physical properties. First of all, we looked at the 
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Figure 18. Position in the galaxy of the three subsaniples of 
clusters. The normal clusters are labeled as filled black points. 
Clusters with flux excess in H and Ka are denoted by red triangles; 
clusters with excess also in I band arc the blue diamonds. 
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Figure 19. Residuals of the chi-square fits in I (bottom), Ks 
(middle) and H (top) versus the estimated ages (left), masses 
(center) and extinctions (right). 



position in the galaxy of clusters affected by a red excess 
(Figure [18)) . They are mainly located in the three starburst 
knots, as one might expect from their young ages. Compar- 
ing Figure [T] and Figure [TS] shows that that many clusters 
with flux excesses are located close to dusty regions. Nev- 
ertheless, we found no clear relation between the measured 
cluster excesses and the corresponding estimated extinctions 
(Figure [191 right column). 

We have looked for further clues to the red excess by 
comparing the deviations from the best SED fits (Am) with 
the derived properties of age, mass and extinction. In Fig- 
ure [19] we plot the measured flux excess in I, H and Ks 
against the recovered ages, masses and extinctions of the 
clusters. We see no correlation between the strength of the 
displacement and any of these properties. The only correla- 
tion with age, visible in Figure [TT] is the absence of clusters 
affected by a red excess with ages older than 40 Myr. This 
result may however be biased by the small number of ob- 
jects in the older age bins. In summary, the lack of any clear 



relations seems to indicate that the red excess could be pro- 
duced from different mecha nisms in diffe r ent clu sters. 

As flrst discussed by iReines et al.l (|2008al ). different 
mechanisms could cause the red excess at such young ages. 
For clusters younger than 5 Myr, a possible explanation 
for the observed I-band excess is the extended red emis- 
sion (ERE) phenomenon. ERE has been observed in a wide 
ran ge of environments, among them star forming regions 
(see lWitt fc Viihl (|2004l ) for a review) and it affects mostly 
optical and near-IR wavelengths, between 7000-9000 A. The 
effect is expected to be strong in presence of strong UV 
radiation fields (as H aro 11 has been shown to possess; 
iGalametz et al.l [2009|) mainly produced by the most mas- 
sive stars in the clusters. However as we have noted, this 
mechanism cannot explain the excess in the infrared, so a 
different physical process needs to be considered. 

Clusters younger than 5 Myr old are still partially em- 
bedded in the original dust cloud. At this stage the com- 
plexity of the cluster environment and the uncertainties in 
the used models are quite large. Hot dust emissions can pro- 
duce a "rise" at IR wavelengths. The strong UV radiation 
field produced by the most massive stars in the cluster is 
expected to be absorbed by dust grains, which eventually 
reemit photons at IR wavelengths depending on the size and 
temperature of the grains. Since Haro 11 has a total FIR 
luminosity that qualifies it as a LI RG, such processes ar e 
undoubtedly working in the galaxy (|Galametz et al.ll2009l ); 
ISmith fc Hancockll2009l). Simi lar explanations were first sug- 
gested bv lReines et al.l ( 2008bj) to explain the red excess they 
observed in the SSCs of SBS 0335-052. New fits to those clus- 
ters (Adamo et al., to be submitted) show that the excesses 
in these clusters can be explained by nebular e mission, a pos- 
sibility which was not properly addressed by iReines et al.l 
|2008b|). On the other hand, at such young stages, it is pos- 
sible that the clusters have two stellar components, one of 
which has been able to shed its dust envelope and contribute 
to the UV-optical emission while the other is still embed- 
ded in the star-forming cloud and still highly extinguished. 
These embedded stars are pre-main sequence objects and 
young stellar objects (YSOs) still surrounded by circumstel- 
lar disks that can produce this IR-bright component. Mod- 
elling is required in order to estimate how massive should be 
the second stellar component in order to produce such large 
displacements in H and Ks band. Moreover, the presence of 
an embedded stellar population contributing only in the IR 
would imply an incorrect mass estimation for the clusters. 

The UV radiation field is expected to soften after a few 
Myr as the most massive stars evolve and finally explode as 
supernovae. This means that in clusters older than 5 Myr 
the I-band excess cannot be explained by the ERE. We also 
expect that for ages > 5 Myr, the embedded phase is over 
and YSOs and pre-main sequence stars have evolved so that 
they do not contribute anymore to the IR fluxes. 

Two further effects could explain the presence of a NIR 
excess in older clusters. (1) Blends of two or more clusters 
may well be present in our data. The pixel size at distance 
of Haro 11 corresponds to ~ 9 pc. If, for example, an older 
cluster is seen in the line of sight towards a very young, 
embedded cluster, the cluster behind would only be visible 
at IR wavelengths, while the estimated age and mass are for 
the UV-optical bright cluster in front. 

(2) Uneven IMF sampling in less massive clusters can 
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also affect the expected NIR flux. Studies of the evolution 
of the stellar isochrones (jMaiz Apellanial2009l . among many 
others) have shown that red supergiants (RSGs) dominate 
the star cluster light in red and near-IR wavelengths as 
early as 6 Myr. ZOl models are constructed by integrat- 
ing the contribution to the cluster's light from all the stars 
populating a well sampled IMF, normalized with respect 
to the mass of the cluster. This method is, however, valid 
only if the stellar population has a total mass of at least 
10^ M0, otherwise uneven sampling introduces important 
changes in the colors of the clusters. In fact, the IMF de- 
scribes a probability distribution of the mass at which stars 
form. This means that, in nature, the IMF sampling is a 
purely stochastic event. For clusters with masses not large 
enough to fully sample the IMF, the colours of the clusters 
are dominated by fluctuations around the expected values 
(see i^ Ccrvifio & Luridiana 2004, Lancon & Mouhcinc 2002, 
iLancon fc Fouesneaull2009l and iFouesneau fc Lanconll2009f ). 
At cluster ages between 6—60 Myr, RSGs dominate the light 
at red and near-IR wavelengths. If the observed cluster hap- 
pens to have a larger number of RSGs than predicted by 
the model, this would be seen as an excess in the observed 
flux at A > 8000 A. Simulated clusters with stochastic IMF 
sampling deviate from the expected evolutionary tracks by 
as much as 1.0 < V ~ H < 5.0 (0.5 < V - I < 2.0), com- 
pared to the expected Q <V - H <1 (-1.0 <V - I < 0.3) 
colours (M. Fouesneau, private communication). 

Finally, a bottom-heavy IMF could in principle explain 
the IR excess over all the age range. However it is diffi- 
cult to imagine why only a fraction of the clusters would 
ma nufacture star s with a different IMF. A recent review 
by iBastian et al.l (|201G| ) , shows observational evidence for 
a universal IMF shape in different galactic and extragalac- 
tic environments. Finally, the total stellar mass necessary to 
produce such an excess would be very high and unphysical 
for "normal" clusters. 



6 GLOBAL STARBURST PROPERTIES AND 
EVOLUTION DERIVED FROM THE STAR 
CLUSTERS 

In Figure [201 we show the position of the clusters in the 
galaxy as function of their estimated ages (colours) and 
masses (circle sizes). The youngest clusters (in blue) are lo- 
cated mainly in close to knot B. The locations suggest a 
starburst propagating from the the older regions C and A 
toward knot B. Older clusters are located throughout the 
main burst regions. The clusters outside the main starburst 
appear to follow a spiral arm-like distribution, but in the 
absence of information on the 3-dimensional distribution we 
do not explore this possibility further. 



6.1 The cluster luminosity function and the 
fraction of missing clusters 

In Figure [TT] we pointed out a deficit of low mass clusters 
in the age range of 1-3 Myr. Using our SSP model, we es- 
timated that clusters at these ages and masses should have 
apparent magnitude in I band in the range 23.0 — 26.0 mag. 
In section 13.21 we have showed that in the crowded, most 
active star-forming regions where we expect to find mainly 
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Figure 20. Spatial distribution of the SCs in Haro 11. The dif- 
ferent circle sizes correspond to four mass ranges (M ^ 10*, 
10* - 10^, 10'^ - 10^ and > WMq). The age ranges are defined 
by colours as indicated in the legend. 



young clusters, detection limits went below 90% already at 
~ 25 mag. Therefore it is possible that a fraction of low mass 
clusters is lost either because they are heavily embedded and 
fall below the detected thresholds imposed by the available 
data. Blending may also affects our sources in that we are 
not able in some cases to resolve single clusters. We consider 
it very unlikely that the missing cluste rs have already bee n 
destroyed (so-called "infant mortalitv". lLada fc Ladall2003l V 
For help with this issue, we decided to investigate the cluster 
luminosity function (CLF). We define the CLF as 



dN{Li) = BL^dLi 



(1) 



where Lt is the luminosity of the clusters in a given wave- 
band i, and B is a constant. In order to determine the slope, 
a, we fit a linear relation in log — log space. 



\ogN{U)^AKh^D, 
with the variable A related to a according to 



A 



-0.4(a-fl). 



(2) 



(3) 



In our specific case we excluded the R filter (F606W) be- 
cause it is heavily contaminated by Ha emission at the 
youngest ages, yielding a shallower CLF. The F814W (I) 
band is instead more suitable for constructing the CLF. As 
we stated in section I3TT1 563 objects were detected in both I 
and R bands. To avoid uncertainties due to underestimates 
of the completeness, we used only clusters detected above 
the 90% completeness limit in I (< 26.5 mag), giving a total 
of 304 objects. In Figure [21] we show the resulting CLF for 
the I band in absolute magnitudes, adopting a distance mod- 
ulus of Af — m = 34.58. The thick black line is the best fit 
to the total distribution, corrected for incompleteness. The 
obtained slope atot = —1.52 ± 0.05, below the typical ob- 
served value, ct ^ —2. Fitting a double power law produces 
small difference between the slope in bright and fainter bins 
so that does not improve the fit. However we noticed in sec- 
tion [3]2] that in crowded regions the 90% completeness limit 
is reached already at ~ 25.0 mag. We therefore performed 
a conservative fit including only objects brighter than 25.0 
mag. The recovered slope is steeper (ai = —1.66 ± 0.03) 
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Figure 21. Observed CLF in tlie F814W filter witli (sofid fine 
liistogram) and withiout (daslied line histogram) corrections for 
incompleteness. The plot shows both absolute and apparent mag- 
nitudes. The thick black line shows the best fit to the complete- 
ness corrected distribution; the slope atot is indicated in the plot. 
Finally, the dashed line in magenta is the best fit line to the dis- 
tribution with lower magnitude limit set to 25.0 mag (the 90 % 
completeness limit for cluster inside the crowded region; see sec- 
tion [3]2]|. The slope ai is labelled in magenta. 
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Figure 22. Cluster formation history (CFH) of the present star- 
burst in Haro 11. The filled black dots connected by the green 
dotted line shows the observed CFH estimated by the total mass 
in clusters more massive than 5 X 10* Mq. The open squares con- 
nected by the dotted red line trace the extrapolated CFH for 
clusters with masses in the range 10^ to 10^ Mq, under the con- 
dition of an ICMF with index —2 and assuming that we know the 
total mass in clusters more massive than 5 X 10* Mq in each age 
bin. 



but still not close to —2. We also tried constructing an 
extinction-corrected CLF, using this time only clusters with 
known extinction. The recovered slope of the distribution, 
a = —1.64 ± 0.02, was still unchanged. 

It is known that cro wding and pr ojection efTects can 
flatten the slope of the LF (|Larsenll2002l ). At the distance of 
Haro 11, and with a pixel resolution of ~ 9 pc in the F814W 
HST filter, some of the clusters are resolved. With an aper- 
ture radius of 4 pixels, our photometry was carried out in a 
region of diameter ~ 72 pc, much larger than typical cluster 
diameters (5 pc). Bright detections may therefore be the re- 
sult of overlapping sources. By substituting, in the equation 
(3), a = —2 and the constant D, given by the assumption 
that the number of objects in the third bin is correct (we 
exclude the two brightest objects, which are knots B and 
C), we made a rough estimation of the missing fraction of 
clusters, using the relation (2). This leads to an estimate of 
85% of objects missing in our analysis. However, this frac- 
tion seems unreasonably high and impossible to reconcile 
with the number of detected clusters and the output of our 
completeness test. While we are certain that blending af- 
fects our data, we cannot give a quantitative estimate of the 
effect on the CLF. Moreover, some of the low mass clus- 
ters with ages between 1 — 3 Myr are still in the embedded 
phase and may fall below the optical detection limits due 
to extinction. Interestingly, in our group's previously pub- 
lis hed CLF study of S Cs in another BCG, ESQ 338-IG04, 
by lOstlin et~ai] (|2003l '). the slope found was -1.8 ± 0.15. 
This galaxy lies close enough (~ 37.4^ Mpc) for blending 
between systems not to be an issue. However, the steepness 
of the CLF is still below the expected value. It is possible 
that the physical environment in BCGs differs from other 
late-type galaxies, perhaps related to their status as merg- 
ing and interacting systems. It has been already noticed that 
a number of dwarf starburst galaxies contain highly lumi- 
nous clusters much brighter that expected from the general 



size o f their clust er population (e.g. NGC 1569. 1 Anders et al 



2004 NGC 1705. lBillett et al.ll2002l: NGC 4 449 



Reines et al 



2008al : SBC 0335-052. iReines et al.ir2008bl ) . iLarsenl (|20oi ) 



refers to NGC 1569 and NGC 1705 as an example of possibly 
biased CLF toward high-mass clusters or in other words that 
CLF in those galaxies may be peculiar. If this is true, then 
BCGs are more similar in this respect to dwarf starburst 
galaxies than to very massive systems. 

6.2 CFR versus SFR in Haro 11 



The st ar formation rate in Haro 11 is very high. lHaves et al.l 
(J2007l ) found a SFR - 19-25Mflyr~ \ using a wide range of 
estima tors. Applying the relation in iRowan- Robinson et al.l 
(|l997h we calculated a total IR luminosity of Lfir = 4.43 x 
10** erg s~^ from the IR luminosity at 60 /xm (taken from 
NED). Using the total IR luminosity of Haro 11 we derived 
the p resent SFR from the Kennicutt's relation (|Kennicuttl 



Il998t ) of SFR = 2I. IM0 yr~\ iii perfec t agreement with the 
estimation made bv lHaves et al.l ((20071). In the present work 
we used a mean value of the SFR = 22 ± ZM^yr'^ . 

Typically, a fraction of the star formation in the galaxy 
occurs in bound clusters, which will possibly survive and 
evolve. The remaining star formation occurs in low-mass 
clusters (analogous to the Milky Way's open clusters), as- 
sociations and more extended systems which would be de- 
stroyed during the first few Myr and disp e rsed i nto the field, 
as predicted by infant mortality. iBastianI (|2008l ) has defined 
the present cluster formation efficiency in a host galaxy as 
r =CFR/SFR, where CFR= Mtot/^t is the present cluster 
formation rate at which the galaxy produces a total clus- 
ter mass Mtot in a giv en age interval Ai. In a recent paper, 
ICoddard et al.l (|2010l ) have found a correlation between the 
value of F and the star formation density Espr of the hosting 
system, confirming that higher SFR leads to the formation 
of more bound clusters. 
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Figure 23. Cluster formation efficiency, F, as function of the 
galactic star formation r ate density, T,sFB,- Th e black diamonds 
are the galaxy sample of lGoddard et alj 1 20101) which were used 
to o btain the best fit po wer-low relation shown by the dashed 
line IIGoddard et al.ll2010l . their equation (3)). At the right hand 
end we show the position of Haro 11 (dot) which fits the relation 
nicely despite its extreme F and SFR values. 



The cluster age distribution for Haro 11 shows that the 
present starburst has lasted for 40 Myr, with a peak of clus- 
ter formation only a few Myr ago. Under the condition that 
a full cluster population forms in 10 Myr, we tried to esti- 
mate how the rate of cluster formation has changed since the 
burst phase began. Taking into account the distance of the 
galaxy and possible blending at the faint end of the CMF, 
we estimated the total mass in clusters more massive than 
5 X IO^'Mq at each 10 Myr age bin, assuming that we detect 
all the clusters younger than 40 Myr with masses above this 
limit. In this way we have constructed the observed clus- 
ter formation history (CFH) in clusters more massive than 
5 X 10* M© (Figure [221). The CFH of massive clusters peaks 
with the youngest cluster population, which has formed in 
the last 10 Myr at a rate of 3.8 MQyr~^. This shows that 
at the present time, r(Msc > 5 x IO^Mq) ~ 0.17, or in 
other words, 17 % of the star formation is happening in mas- 
sive clusters. The CFR declines toward lower rates at larger 
lookback times and reaches a value of 0.08MQyr~^ when the 
cluster population formed around 40 Myr ago. This decline 
may also indicate lower star formation rates at this epoch. 
However the massive cluster formation in the galaxy during 
the whole period has been extremely efficient. 

In order to constrain a total cluster formation rate, we 
estimated the missing fraction of mass in low mass clusters 
(10^ < Mqc < 5 X IO'^A'/q), taking an initial cluster mass 
function (ICMF) with slope —2 and assuming that the total 
cluster mass in the range 5 x 10* — IQ^ Mq is known. The 
extracted total CFH is also shown in Figure 1221 The total 
CFR at ages below 10 Myr is 8.31M0yr~^, which means 
cluster formation efficiency F ^ 0.38 ± 0.6 or ~ 38 %. This 
fraction is extremely high compared to the value observed 
in cluster systems in other galaxies (see Table |4| , and al- 
most a factor of 10 hig her than in the SMC (F = 2 - 4%; 
Gieles fc BastiaJl 2008l ) and the solar neighborhood (~ 5%; 
Lada fc L ada 2003). We reach the same result if we esti- 
mate the present CFR from the most massive cluster, Mmax, 



in the age bin 1-10 Myr. In fact ICieles fc BastianI (|2008l ') 
showed that Aftot — lOAfmaa;, for a constant cluster forma- 
tion rate and a power-law CIMF with index —2, and assum- 
ing that the maximum cluster mass formed by the galaxy, 
Mup is bigger than the Mmax observed in the system. If we 
assume that the CFR is constant over the last 10 Myr and 
that the most massive cluster formed^j is knot B with a 
mass of 8.35 x lO^'M©, the present CFR is 8.35Af0yr"\ 
very similar to the value estimated above. 

If blending affects our cluster mass estimates then we 
may be overestimating t he CF R. Applying the relation 
found by ICoddard etall (|2010h between r(%) and Esfr, 
we can estimate the expected value of F for the SFR ob- 
served in Haro 11. The starburst in this galaxy covers an 
area of ~2.34 kpc^ so that the Ssfr = 1.28M q yr"^ kpc"^. 
From equation 3 in ICoddard et al.l (|20ld ) we find that 
r = 31.0 ± 5% which is in quite good agreement with the 
previously fo und value. In F i gure 1 231 we show the relation 
presented bv lCoddard et al.l (|2010r ) and the galaxy sample 
they used to derive it, with the addition of the position of 
Haro 11. We observed that in spite of the extreme cluster 
formation environmen t in Haro 11, the correlation found by 
ICoddard et all (|2010l ') between F and density of SFR is still 
valid, supporting the theory of a close connection between 
the global properties of the host galaxy and its star cluster 
population. 

In Table|4]we summarise the properties of the Goddard 
et al. sample together with Haro 11. For some of the objects 
the authors gave the total mass contained in clusters younger 
than 10 Myr. We used these values to estimate at which 
fraction of the total mass of the system they do correspond, 
as indicated in the last column of Table 3] For Haro 11 we 
used only the total observed mass in clusters more massive 
than 5 x 10*^/0. The comparison between the values for 
different galaxies nevertheless gives an idea of the strength 
of the cluster formation process in Haro 11. 



6.3 Discussion 



In a recent work. lMelena et al.l ((20091) studied the properties 
of clumpy star-forming regions associated with bound clus- 
ters in nearby dwarf starburst galaxies. They found that, 
independently of their environmental conditions, the for- 
mation of those clumps is given by random sampling of 
a mass function of the form dN/dM oc M~^ . We can 
think the mass function as a statistical relation which is 
stochastically realized by star formation in clusters. The 
galactic environment limits the range of the stochastic ac- 
tion. For example, the galactic environmental conditions can 
limit the maximum possible m ass at which a cluster can 
be formed (|Gieles et al.ll2006al ). The galactic environment 
may also lead to a characteristic mass, Mc, above which 
a CIMF shaped like a Schcchtcr function decline s expo- 
nentia lly (|Larscn 2009; Gjclcs 200y). In particular, iLarserJ 
(|2009h noticed that extreme star-forming environments show 
no characteristic truncation at 10^ < Mc < 10^ because 
such galaxies are actually able to form very massive clusters 



^^ The very massive region C is excluded due to the uncertainties 
in its SED fit. 
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Table 4. The cluster population of Haro 11 compared to the galaxy sample of iGoddard et alj l|201Clh . We show the SFR, r(%) and 
total mass in clusters younger than 10 Myr and more massive than lO^M© (NGC 1569, M83, NGC 6946) or 5 X 10* Mq (NGC 3256 
and Haro 11), MysCi the total hosting galaxy mass, Mg^i, and the fraction of the total cluster mass with respect to t he galactic mass, 
/ =MYSc/M j,ai. The dat a of t he galaxies NGC 1569 , M8 3, NGC 6946, and N GC 3256 are fromlCoddard et al.1 | |2010|'I unless otherwise 
specified, (a'l flsraell digSSD : (b) iLundgren et al.1 ||2004| '): (c) IWalsh erZI | |2002|'I : (d'l lFeast fc RobertsonI l ll978l) . 



galaxy 



SFR(Mo yr-i) r(%) MyscCMq) Mg,i(M0) 



/ 



NGC 1569 


0.3626 


13.9±0.8 


3.52x10^ 


4.6x10** (a) 


0.0008 


M83 


0.3867 


26.7±t:g 


7.24x10^ 


6.1xl0i0{b) 


1.5x10-5 


NGC 6946 


0.1725 


12.5±i| 


1.284x10^ 


5x101''(c) 


2.6xl0-« 


NGC 3256 


46.17 


22.9±^i 


1.66x10'' 


- 5 X lOiO(d) 


0.0003 


Haro 11 


22 ±3 


38±6 


3.82x10'^ 


~ I X 10^0 


0.004 



(M > lO^V lLarsenl (|2009|) presented two hypotheses regard- 
ing the physical mechanism responsible for the formation of 
massive clusters. In the first, such clusters are formed by 
super-GMCs, with masses of at least 10^ — IQ^Mq. The sec- 
ond requires a higher star formation efficiency during the 
collapse of the GMC, which in turn implies compression of 
the gas to higher densities. In the Milky Way it has been 
observed that GMCs have low de nsities and form clu sters 
with low star formation efficiency (JLada fc Ladall2003l l pro- 
ducing low mass clusters. In quiescent massive regular spi- 
rals like the Milky Way the shear due to the spiral patterns 
would then favours the fragmentation of GMCs, preventing 
very massive cl usters from forming . On the other hand, as 
pointed out bv lBillett et al.l ()2002l ). in dwarf starburst sys- 
tems the gravitational instabilities acting on the GMCs are 
instead much stronger than the shear from the irregular in- 
tergalactic medium, allowing the observed massive clusters 
to form. 

In Haro 11 we have found many massive clusters, over 
35% of which have masses above IQ^Mq. In the previous sec- 
tion we showed that the efficiency with which the galaxy has 
formed clusters, in particular massive clusters, is higher than 
in merging systems like NGC 3256. Haro 11 has an irregular 
morphology and is less massive than normal spiral galaxies. 
It seems likely that a previous merger with a gas-rich system 
has favoured the inflow and compression of gas, triggering 
the formation of the massive clusters. In this way, the galaxy, 
despite its peculiarities, easily fits into the wider picture of 
a sequence in star formation intensity stretching from quies- 
cent to very active modes according to Goddard's relation. 
The age distribution of the clusters confirms th e present es- 
timate d starburst age of ~ 35 Myr found by lOstlin et al.l 
l|200ll ). The starburst has propagated rapidly through the 
entire galaxy as showed in Figure [20l consuming the avail- 
able fuel and is possibly now settling into a period of more 
modest activi ty, as attested by th e low H I mass content 
in the galaxy iBergvall et al.l l200d . From the available H I 
mass, the remaining life time of the present burst phase is 
~ 5 Myr. 

A final important issue is raised by the similar- 
ity b etween Haro 11 an d LBGs observed at high red- 
shift ijOverzier et al.ll2008h . Nearby Lyman break analogues 
(LBAs) are attractive objects for studying star formation 
modes a nd morphology i n det ail which is impossible at high 
redshift. lOverzier et al.l (|2009l ) have for example shown that 
star formation in LBAs is both very active and concentrated 
in clumpy regions which may be composed of supermassive 



star clusters. In Haro 11, we found a CFR of 38%, somewhat 
higher than locally observed. This result suggests that the 
CFRs in LBGs could be almost a factor of 10 higher than 
in the Local Group, implying that at high redshift galaxies 
formed more numerous and more massive clusters than we 
observe in the local Universe. However investigation of other 
systems similar to Haro 11 are needed before we can draw 
firm conclusions. 



7 CONCLUSIONS 

Thanks to the deep HST data stretching from UV to IR 
wavebands, we have carried out a photometric study of the 
clusters in the starburst regions of Haro 11. We have con- 
strained the physical properties of 185 clusters through the 
analysis of their integrated SEDs. We summarize the main 
results here. 

(i) More than 60 % of the star clusters in Haro 11 are 
younger than 10 Myr old. The present starburst phase 
started less than 40 Myr ago with a peak as recently as 3.5 
Myr. Only 16 clusters have estimated ages greater than 40 
Myr and, of these, half are older than > 100 Myr. Haro 11 
therefore provides a unique opportunity to investigate the 
cluster formation and evolution in a starburst system in its 
earliest phases. For ages between 1 and 3.5 Myr the extinc- 
tion range is quite wide, consistent with a complex cluster 
formation process and the gradual dissolution of the parent 
clouds. 

(ii) The clusters are quite massive, with masses between 
10^ - 10® Mq. There are around 60 estimated SSCs, 7 of 
them with masses larger than l(f Mq, and two very massive 
knots (B and C) with estimated masses ^ lO' Mq which we 
could not resolve into smaller clusters. Very young clusters 
(1-3 Myr) with masses around 10'' — lO'^M© are absent from 
our data even though our detection limits should have al- 
lowed their detection. These missing clusters could still be 
embedded in their parental clouds, invisible in optical wave- 
bands, or blended with other clusters in the most crowded 
regions. 

(iii) We have discovered that a fraction of the clusters 
exhibits a fiux excess compared to plausible SED models at 
wavelengths > 8000 A. We have noted that fitting those clus- 
ter SEDs affected by the red excess with normal SSP models 
leads to overestimation of the clusters ages and masses. In 
several cases, we have obtained acceptable fits only at bands 
blueward of I, highlighting the potential impact of the phe- 
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nomenon on estimations of the cluster properties. Of the 
clusters affected by the red excess, none are older than 40 
Myr. 

(iv) We have discussed possible causes for the red excess. 
If the age of the cluster is the key property, then the red 
excess in clusters younger than 5 Myr could be due to the si- 
multaneous presence of very young UV-bright massive stars 
and the dust left from the GMCs. The ERE mechanism 
could then originate the excess at I waveband while hot dust 
emission produces a "rise" at IR wavelengths, depending on 
the size and temperature of the dust grains. Another possible 
mechanism at such young stages could be the presence in the 
cluster of two stellar components, one of which has shed its 
dust envelope, contributing mainly in the UV-optical wave- 
lengths, while the other component is still deeply embedded 
in its star-forming cloud and highly extinguished. Pre-main 
sequence objects and young stellar objects (YSOs) still sur- 
rounded by circumstellar disks can produce this IR-bright 
component. Further modelling will be required in order to 
estimate how massive the second stellar component needs 
to be in order to produce such large displacements in the 
H and Ks bands. The presence of such embedded stellar 
populations emitting only in the IR would also imply that 
the mass of these clusters will be underestimated. At ages 
> 5Myr we expect that the embedded phase is over; we 
have suggested two further mechanisms to explain the NIR 
excess for ages older than this. First, evolved clusters may 
be observed in the line of sight towards younger, embedded 
ones. In such cases, the cluster behind would only be visi- 
ble at IR-bands, while the age and mass estimated form UV 
and optical bands correspond to the cluster in front. Second, 
red supergiants, which dominate the light at red and near- 
IR wavelengths, could "redden" the observed colours of the 
clusters at ages as early as 6 Myr. 

(v) For the fir st time, after an early attempt by 
lOstlin et al.l (|2003r ). we have analysed the cluster luminosity 
function in a BCG. Haro ll's CLF is shallower than ob- 
served in other galaxies. Contamination by blending could 
have this effect on the CLF, but we didn't find any evi- 
dence of overestimation of the CFR, using the more massive 
SSCs. Obviously, at the distance of the galaxy and with the 
present available data, we cannot avoid crowding and pro- 
jection contaminations. 

(vi) More than the 30% of the star formation in Haro 11 
has happened in bound clusters which can survive and evolve 
along with the galaxy. This fraction is higher than observed 
in the Milky Way and other merging galaxies like NGC 3256, 
leading to the conclusion that the environmental conditions 
in this BCG have favoured the formation of massive super 
star clusters. Due to the similarities between Haro 11 and 
the high-redshift LBGs, we suggest that cluster formation 
rates at higher redshift may also be enhanced. 
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APPENDIX A: CHECK ON THE 
CALIBRATION OF THE F160W FRAME 

We performed 2 different tests on the HST IR-frame in order 
to verify the goodness of the calibration. First, we compared 
photometry output of few bright sources between the final 
reduced frame and one single raw frame. At the same aper- 
ture radius, sources on the raw frame showed to be ~ 0.15 
mag brighter than the final science frame. However this scat- 
ter can be explained by sky background and bias effects cor- 
rections performed on the final frame. In the second test, we 
used an available HST archive Af/CS/FlOOW exposure of the 
MRK 930 target. Because of the analogy between the NIC3 
data set of Haro 11 and MRK 930 we decided to compare the 
photometry outputs of MRK 930 NIC2 and NIC3 F160W 
frames. The NIC2 frames were reduced with MULTIDRIZZLE. 
For small aperture radii, the photometry of sources in the 
NIC2 frame were brighter than the NIC3 ones with a similar 
constant scatter. At larger radii, the scatter became rapidly 
negligible. As expected this behavior is due to the better 
pixel resolution scale and PSF sampling of the NIC2 cam- 
era. However because the photometry converged at the same 
values this shows the correctness of the calibration, showing 
at the same time the importance of an accurate aperture 
correction to infinity. 
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Table 1. Photometric properties of the 4 clusters showed in Figure [4] Magnitudes are in Vega system throughput. 



ID 



F606W 



F814W 



F140LP 



F220W 



F330W 



F435W 



F550M 



F160W 



Ks 



6 22.97±0.1 

43 21.49±0.11 
87 22.32±0.19 



22.67±0.07 
20.81±0.05 
22.16±0.17 



113 22.49± 0.1 22.09±0.05 



20.75±0.06 
21.70±0.11 
19.84±0.12 



21.25±0.07 
21.73±0.13 
20.51±0.14 



22.10±0.08 
21.59±0.09 



23.12±0.05 23.21±0.06 21.24±0.42 20.18±0.66 

22.54±0.08 22.33±0.07 19.07±0.4 17.71±0.64 

22.60±0.15 22.42±0.13 - 19.21±0.66 

23.29±0.06 22.85± 0.06 21.11±0.42 20.12±0.66 
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by the author. 
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